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ABSTRACT 

We find that the formation of MWC 656 (the first Be binary containing a black 
hole) involves a common envelope phase and a supernova explosion. This result sup¬ 
ports the idea that a rapidly rotating Be star can emerge out of a common envelope 
phase, which is very intriguing because this evolutionary stage is thought to be too 
fast to lead to significant accretion and spin up of the B star. We predict ~ 10-100 
of B-BH binaries to currently reside in the Galactic disk, among which around 1/3 
contain a Be star, but there is only a small chance to observe a system with parame¬ 
ters resembling MWC 656. If MWC 656 is representative of intrinsic Galactic Be-BH 
binary population, it may indicate that standard evolutionary theory needs to be re¬ 
vised. This would pose another evolutionary problem in understanding BH binaries, 
with BH X-ray Novae formation issue being the prime example. The future evolution 
of MWC 656 with a ^ 5Mq black hole and with a ^ 13Mq main sequence companion 
on a ^ 60 day orbit may lead to the formation of a coalescing BH-NS system. The 
estimated Advanced LIGO/Virgo detection rate of such systems is up to 0.2 yr“^. 
This empirical estimate is a lower limit as it is obtained with only one particular evo¬ 
lutionary scenario, the MWC 656 binary. This is only a third such estimate available 
(after Cyg X-1 and Cyg X-3), and it lends additional support to the existence of so 
far undetected BH-NS binaries. 
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1 INTRODUCTION 

We know at present 184 X-ray binaries consisting of a B e 
star and a compact object - Be XRBs JZioIkowskil l2014h . 
Until the previous year, whenever the nature of the com¬ 
pact component was determined (in 119 systems), it was 
always a neutron star. Not a single Be system containing 
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a black hole was found during 40 years of observations of 
Be XRBs. Last year th e first such system was finally found 
bv ICasares et al.1 (l2014l ). This discovery motivated us to in¬ 
vestigate the possible evolutionary scenarios leading to the 
formation of similar systems. 

1.1 Be stars 

Be stars are massive, not substantially evolved, main se- 
quen ce stars of spectral types BO-AO with Balmer emission 
lines JPorter fc Riviniual2003h . This range of spectral types 
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corresponds roughly to a mass range of about 3 to 18 Mq. 
The emission lines (which give the name to this class of 
stars) originate in an outflowing viscous disc (excretion disc) 
around the star. Snch a disc is very similar to the well known 
viscous accretion discs, except for the changed sign of the 
rate of the mass flow. The excretion discs evolve dynamically 
on a time scale of a few years to few decades. In the course of 
this evolution, th e disc unde rgoes a global one-armed oscil¬ 
lation instability (lKatolll983fl . manifesting itself in the form 
of the well observed, so called, V/R variability. This insta¬ 
bility (progressing density waves) leads eventually to the 
disruption of the disc. A disc-less phase (with no emission 
lines) follows the n, until the disc refills ag ain (which takes 
years to decades (IPorter fc Riviniua 2003lll. The excretion 


discs were suc c essfully modeled jHummel fc Vranckerj 


Okaza^ 19961: Hummel fc Hanuschik 19971: Okazaki 


1995 


1997 


Porter! 19991 : iNegueruela fc OkazaS " 2000l . l2001^ arid these 

modelings helped to explain many observed properties of 
Be stars (both solitary ones and those that are members of 
binaries, in particular, X-ray binaries). 


1.2 Fraction of Be stars among B stars 

At first, we would like to comment on the notion that the 
same star might show a Be phenomenon (and so be a Be 
star) over some intervals of time and not to show it (and so 
be a “normal” B star) over some other intervals. It is quite 
likely, that a “normal” B star in the course of its evolution is 
spun up and develops a Be phenomenon. However, it is clear 
that there are B stars (majority of them) that are not Be 
stars and that there are Be stars that show Be phenomenon 
during all time we observed the m. They exhibit disc-les s 
phases, which may last decades ijPorter fc Riviniua l2003ll . 
Such phases are part of a Be phenomenon and the stars 
do not stop being Be stars (at least, this is the accepted 
convention). 

After this clarification, let u s estimate h ow large a frac¬ 
tion of all B stars are Be stars. lAbtl (Il987l i found that Be 
stars comprise 18% of the B0-B7 stars in a volume-limited 
sample of field stars (with a maximum Be fraction for the 
spectral types B3-B 4). A similar result was obtained by 
IZorec fc Briod (Il997l l who found the fraction of Be stars to 
be 17% for the Galactic field stars (this was a mean value; a 
maximum of 34% was found for the spectral type Bl). The 
percentage of Be st ars was estimated also for stellar clusters. 

I Keller et al.l lll999ll investigated the frequency of Be stars in 
six young clusters in the Magellani c Clou ds. They found 
a range of 13 to 34%. iMaeder et al.l (Il999li investigated 21 
clusters in the interior of the Galaxy, the exterior of the 
Galaxy, the LMC and the SMC and found the fractions of, 
respectively, 1 1, 19, 23 and 39%. The s imilar investigation 
carried out bv lWisniewski fc BiorkmanI (l2006l l brought sim¬ 
ilar results (the ranges of Be fractions were 9 to 39% for ear¬ 
lier typ e (B0-B3) Be sta rs and 3 to 32% for later (B4-B5) 
types). iMcSwain fc Giei ll2005! l analyzed 48 open clusters 
and found the mean (for all clusters) percentage of Be stars 
equal to 7.1% (with a maximum of ~ 11% for spectral types 
B2-B3). They compare it with Abt’s value and attribute 
the difference to the selection effects in Abt’s estimate (a 
bias t owards earlier B spectral tvpesl. [Pabregat fc TorreionI 
(l200(tl analyzed seven “Be rich” clusters in Milky Way and 
in the Magellanic Clouds. They found very high percentage: 


21 to ^ 50%. Anothe r “Be r ich” cluster was investigated by 
iMarco fc Negueruelal 1 I 2 OI 3 !) who found very high (~ 40%) 
fraction of Be stars close to the turn-off (spectral type Bl) 
but very few Be stars for later spectral types. 

To summarize, the fraction of Be stars among B stars is 
about 20-30%, generally increasing for earlier spectral types. 
If we consider Be X-ray binaries, one should remember that 
Be stars in X-ray binaries have somewhat earlier spectral 
types than solitary Be stars (see the further text), which may 
indicate that the factor /se (the fraction of X-ray binary 
containing a Be star) should be somewhat higher - perhaps 
30%. 


1.3 The origin of fast rotation of Be stars 


There is little doubt that all distinct properties of Be stars 
are related to the presence of an outflowing excretion disc. 
There is also little doubt that the presence of thes e discs 
is related to the fast rotation of these stars. IStruvel (Il93lfl 
suggested that this rotation is very close to the critical 
(or break-up) e quatorial velocity. Later, a canoni cal view 
was established IPorterll 19961 : iGhauville et al.ll200ll l accord¬ 
ing to which the rotation is significantly subcritical with 
equatorial velocity equal only 70-80% of the critica l ve- 
locity. Howeve r , more recently iTownsend et al.l ll2004ll and 
lEkstrom et al.l ll2008ah gave arguments indicating that the 
rotation might be indeed very close to critical, with equato¬ 
rial velocity smaller only by a few percent (and not 20-30%) 
than the critical one. Such fast velocity makes the forma- 
tion and maintaining of the excretion disc much more likely 
llGranada et al.ll201lj '). 

As for the origin of fast rotation. iMartavan et ah] (l2006fl 
indicated that Be stars are born with higher initial (on 
ZAMS) rotation than other B stars. It seems, however, that 
they are not born as Be stars from the very beginning. 
Rather, the higher initial rotation facilitates the action of 
mechanisms that later spin up these stars to nearly critical 
rotation. Two such major mechanisms were considered. 

One of them is the evolutionary spin up during the 
Main Sequence (MS) evolution. The reason for the spin 
up is the significant decreasing of the moment of inertia 
of the star during this phase of evolution. This explana¬ 
tion might be supported by the fact the Be phenomenon 
seems to be ass ociated with the second half of the MS evolu - 
tion of B stars llMcSwain fc Gie^l2005l : [Premat et al.ll200^ 1. 
The first me chanism was discussed and modeled by differ¬ 
ent authors llMevnet fc Maederll2005l : lEkstrom et ai.l[2008al : 
iGranada et al.ll2013fl . The general conclusion of these mod¬ 
elings is that evolutionary spin up is sufficient to explain the 
Be phenomenon. 

The second mechanism is the spin up due to ac¬ 
cretion in a binary system. Initially this scenario was 
proposed for the formation of Be -X-ray binaries (e.g., 
iRappaport fc van den Heuvell lll982ll '). This scenario can 
also account for single Be stars, due to disrupted binary 
systems or binary mergers. First esti mates of the impor ¬ 
tance of th i s scenario were provided by Waters et al.l lll989ll ; 

I Pols et ^ (I 1991 I I ; jPortegies ZwartI (Il995l l . More recent sim¬ 
ulations accounting for the actual spin u p pro cess and merg¬ 
ers were performed by Ide Mink et al.l (l2013l l . The general 
view of the advocates of the binary mechanism is that this 
scenario can account for the majority of the Be stars. For 
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example iMcSwain fc Gie3 (l2005l l concluded that more than 
70% of all Be stars had to be spun up in the p roce ss of accre¬ 
tion in binary systems. Ide Mink et alJ (1201.11 1 an d IShao fcTil 
(|20141 both conclude that all Be stars can be accounted for 
by the binary evolution. Especially, if mergers of short pe¬ 
riod contact binaries are taken into account. 

The summary of this controversy is difficult. It seems 
possible that both mechanisms are at work and that, at 
present, it is not possible to estimate reliably the relative 
importance of each formation channel. 


1.4 Be stars in X-ray binaries 

Be stars in X-ray binaries are, in many respects, similar to 
isolated Be stars. They have excretion discs which develop 
one-armed oscillations and display V/R variability. The dy¬ 
namical evolution of their discs includes the disruption of 
the disc and the following disc-less phase. 

However, there are also notable differences. The first 
concerns the spectral types. Be stars in X-ray binaries have, 
on average, earlie r spectral types than isolated Be stars 
dNegueruelal 1 1998 1 . The range is 09 to B3 as opposed to 
BO to AO for isolated Be stars. It indicates that also masses 
of Be stars in X-ray binaries are somewhat higher (perhaps 
6 to 24 Mq instead of 3 to 18 Mq). The second major differ¬ 
ence concerns the excretion discs. The discs in X-ray binaries 
interact with the compact companions (almost always neu¬ 
tron stars). This interaction leads to the appearance and the 
growth in the certain locations in the disc of the resonances 
between Keplerian frequencies of the disc matter and the or¬ 
bital frequency of the neutron s tar. These resonances lead to 
the tidal truncation of the dis c dArtvmowicz fc Lubow|[l9^: 
Neeueruela &: Okazakil 120001 . 1200 ll : lOkazaki fc Neeueruelal 
2001^ . Tidal truncation makes the discs in Be X-ray binaries 
smaller and denser than the discs around isol ated Be stars 
dReig et al. I I 1997 I : iNegueruela fc OkazakilliooH l. 


1.5 Be X-ray binaries 

Be X-ray binaries (Be XRBs) are the most numerous class 
among high mass X-ray binaries. We know at present 184 
Be XRBs and only about 60 other high mass X-ray binaries 
(both Gal actic and extragal actic systems are included in this 
statistics: IZiolkowskil d2014l ')l. In 119 of Be XRB systems the 
X-ray pulsations are observed, confirming that the compact 
component must be a neutron sta r. The pulse periods are 
in the range of 34 ms to 1400 s dZiolkowski fc Belczvnskil 
d201llH . Until the previous year not a single Be system con¬ 
taining a black hole was found. The Be XRBs are rather 
wide systems (orbital p eriods in the range of ~ 10 — 1180 
days; IZiolkow^ d2014l H. The orbits are frequently eccen¬ 
tric. A compact component accretes from the excretion disc 
of a Be star (earlier known as the equatorial wind of a Be 
star). 

The X-ray emission from Be XRBs (with a few excep¬ 
tions) is of a distinctly transient nature with rather short 
(days to weeks) active phases separated by much longer 
(months to tens of years) quiescent intervals (a typical flar¬ 
ing behavior). There are two types of flares, which are clas¬ 
sified as Type I outbursts (smaller and roughly regularly re¬ 
peating) and Type II outbursts (larger and irregular). This 


classification was first defined by IStella et al.l dl986ll . Type 
I bursts are observed in systems with highly eccentric or¬ 
bits. They occur close to periastron passages of a neutron 
star. They are repeating at intervals Porb. Type II bursts 
may occur at any orbital phase. They are correlated with 
the disruption of the excretion disc around Be star. They 
repeat on time scale of the dynamical evolution of the excre¬ 
tion disc (~ few years to few tens of years). This recurrence 
time scale is generally m uch longer than the orbital period 
dNegueruela et aDl200ll ~l. 

Be XRBs systems are known to contain two discs: an ex¬ 
cretion disc around the Be star and an accretion disc around 
the neutron star. Both discs are temporary: the excretion 
disc disperses and refills on time scales ^ few years to few 
decad es (dynamical evolution of the disc llPorter fc Riviniua 
l2003lH . while the accretion disc disperses and refills on time 
scales weeks to years (which is related either to the or¬ 
bital motion of a neutron star on an eccentric orbit or to the 
disruption episodes of the excret ion disc). Formation of the 
accre tion discs was ana lyzed by (iHavasaki fc Okazakill2006h 
and llCheng et al.ll20l3 l. 

The more detailed descri ption of the pr o pertie s 
of Be XRBs is given, e.g. 


Ziolkowski fc Belczvnski 


ption of the pr o pertie s 
jNe guemela^^tall (2001); 


Ziolkowskl__([2b(^h^_ Belczroski fc Ziolkowskil 1 2009 1: 


( 2011^ : [ileigl 1 2oT3) and refer¬ 
ences therein. 

The fact that we observe over one hundred neutron 
star Be XRBs and not a single black hole Be XRB, became 
known as the problem of the missing black hole Be XRBs. 
Trying to explain the reasons for which we do not ob serve 
black hole Be XRBs, iBelczvnski fc Ziolkowskil (l2009fl car¬ 
ried out stellar population synthesis calculations aimed at 
estimating the ratio of neutron star to black hole Be XRBs, 
expected on the basis of the stellar evolution theory. The 
results of their calculations predict that for our Galaxy the 
expected ratio of Be X-ray binaries with neutron stars to the 
ones with black holes Pns/bh should be, most likely, equal ~ 
54. Since we know 48 neutron star Be systems in the Galaxy, 
then it comes out that the expected number of black hole 
systems sh ould be just one. It seems that this system was 
just found dCasares et al.ll20141 . 


2 MODELING 

2.1 The StarTrack code 

We use the StarTr a ck p opulation synthesis code 
l|Belczvnski et all I2OO2I . l2008lf to generate a popula¬ 
tion of binaries in the Galaxy . The co de is based on revised 
formulas from iHurlev et al.l (l2000ll . among other with 
updated wind mass loss prescriptions, and calibrated tidal 
interactions, physical estimation of donor’s binding energy 
(A) and convection driven, neutrino enhanced supernova 
engines. A full description of the code can be found in the 
papers mentioned above. Here we concentrate only on these 
aspects that are important from the viewpoint of modeling 
used in this study. 

The initial parameters of systems simulated with the 
StarTrack are randomly chosen from the following distribu¬ 
tions: 

- three component broken power-law initial mass function 
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(IMF) with slope of —1.3 for initial mass Mzams = 0.08- 
0.5Mq, —2.2 for Mzams = O.S-lM pi, —2.7 for Mzams = 1- 
150M© (iKroupa fc Weidn^l2003l l 

- flat mass ratio d istribution in range q = 0 — 1 

l|Kobulnickv et ^l200fil 

- flat i n the loga rithmic distribution of initial binary sep¬ 
arations (lAbtlll98al oc 1/a in range from a minimum value, 
to prevent stars from filling their Roche lobes at zero-age, 
up to 1 O®J ?0 

- thermal-equilibrium distribution of eccentricities 
dPuauennov fc Mavorj|l99lh H(e) = 2e in a range 0-1. 

The binary fraction is assumed to be 50% (i.e., 2/3 of stars 
are in binaries). We no te that recently m easured initial dis¬ 
tributions for O stars dSana et al.ll20lj l are different from 
those employed in our study. However, it was demonstrated 
that a change of distributions from the ones used here to the 
new ones does not signihcantly affect predictions for progen¬ 
itors of double compact objects (de Mink & Belczynski 2015, 
in prep). 

Among the most important physical mechanisms driv¬ 
ing binary evolution is the common envelope (CE) phase. 
The CE is very efficient in creating close binaries. The out¬ 
come of this evolutiona ry phase can b e described by the 
energy balance formula (IWebbinkll 1984h : 

/ GMdon.f Afacc GMdon.iAfacc \ GMdon.i Afdon.env 

V 2Af ~ 2Ai J ^ ARdon.iob 

( 1 ) 

where Mdon and Mace are masses of donor and accretor, 
respectively; Mdon.env is mass of the envelope of the accre¬ 
tor; A is a binary separation; i?don,iob is a Roche lobe ra¬ 
dius of the donor at the beginning of mass transfer. Index 
i/f indicates the initial/final (before/after CE) value of a 
given quantity. The A parameter describes the binding en¬ 
ergy of the envelope of the dono r and in the current version 
of code we use the “Nanjing” A dXu fc Lil (l2010bl '). IXu fc Lil 
((2 oTo 3)) with sp ecific implementation into the StarTrack 
code described in iDominik et aO (120121 '). The acs parame¬ 
ter describes the efficiency of the transfer of orbital energy 
into the envelope. We allow for large variation of this pa¬ 
rameter. 

Another important phase in the past evolution of bina¬ 
ries is the core-collapse/supernova (SN) explosion. Due to 
its potential asymmetry, the new born compact object (NS 
or BH) may receive a nat al kick. According to the observed 
velocities of radio pulsar dHobbs et al.ll20(^ l. we choose the 
maximum kick velocity from the single Maxwellian distri¬ 
bution with a = 265 km s“^. This value can be reduced 
depending on the fallback factor (/fb), which describes per¬ 
centage amount of matter ejected during SN and accreted 
back onto the compact object: 

14 = I4nax(l — /fb) (2) 

This prescription is used for neutron stars (NS) and black 
holes (BH), but most of the former receive full natal kicks, 
with exception of electron capture SN (ECSN) for which we 
adopt no natal kicks at all. The full description of double 
co mpact objec t form ation and rationale behind i t is given 
bv iFrver et al.l d2012l ') and iBelczvnski et al.l d2012l '). 


2.2 The standard model 

In the standard model we employ energy balance for CE evo¬ 
lution with fully efficient transfer of orbital energy to the en¬ 
velope energy (a = 1). The maximum natal kicks velocities 
are drawn from the Maxwellian distribution with a — 256 
km s“^. Fallback factor ffi, varies in range 0-1 (from full kick 
for 0 to no kick for 1). The SN explosion mechanism for the 
standar d model is a con vection driven, neutrino enhanced 
engine (iFrver et al.ll20l3 j. The explosion occurs within the 
first 0.1 — 0.2 s (so called ’’rapid” explosion). This engine re- 


produces (IBelczvnski et al. 20121 

the observed Galactic X 

rav binary mass gap (Ozel et al 

I 2 OIOI: iBailvn et al.lll998ll 

All results that we present are 

obtained for the metallic- 


ity typical of the Solar neighbourhood Z = 0.02 and pre¬ 
sented for a specific assumption on CE outcome: all donors 
beyond main sequence are a llowed to survive CE (but see 
IBelczvnski et al] (l2007l . l2010lj for an alternative scenario for 
Hertzsprung Gap donors). Eull des cription of the standa rd 
evolutionary model can be found in iDominik et all (l2012|j . 

We only evolved binaries with primaries in mass range 
6 -I 5 OM 0 and secondaries in mass range l-lhOM©, as lower 
mass binaries may have only a very little chance to produce 
binaries containing a BH. We only follow stars in the Galac¬ 
tic disk for which we assume a constant star formation at 
the level of 3.5 M© yr“^ over 10 Gyr^. We have evolved 10® 
binaries which gives us one full realization of the Galactic 
disk. We then use Monte Carlo techniques to evaluate the 
probability of catching each system during B-BH phase at 
the current time in the Galactic disk. A fraction /bs of B-BH 
binaries belong to the interesting group of Be-BH systems 
(see Sec. 11.211 . 


2.3 Variations on the standard model 

We also check the sensitivity of our predictions to parame¬ 
ters that are important for the formation and evolution of 
B-BH binaries. 

In models Vi and 14 we change the acs value to 0.1 
and 5 respectively. In model V 3 all BH do not receive any 
natal kicks during SN explosion. In model 14 we change 
the mechanism of SN explosion from the “rapid” explosion 
to the “delayed” one. In delayed SN engine, the explosion 
can occur as late as Is after bounce. This scenario, unlike 
the “rapid” one, produces a continuous mass spectrum of 
compact objects, and the observed mass g ap must then be 
a res ult of some observational bias (e.g., iKreidberg et al.l 

I2012D . 


^ iDiehl et al.l _ 1 20061^ . iMisiriotis et alj ll2006h . 

iRobitaille Sz Whitney ll2010h find a total star form ation rate 
to be 4, 2.7, ~ 1 Mq yr“^, respectively. A study by iKennicuttI 
lll998h suggests SFR value in the range 0-10 Mq yr“^. Ad¬ 
ditionally, if we assume a constant SFR over the period of 10 
Gyr then the Galactic disc mass provided by iMcMill^ ll201lh 
divided by this period gives an average SFR 6.5 Mq yr“^. 
Thus, the value of SFR chosen for our study (3.5 Mq yr~^) is 
within reasonable limits. 
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2.4 Definition of B-BH and MWC 656 - like 
systems 

Due to large uncertainties in rotation physics we choose not 
to use rotation as a qualiher for the Be phenomenon. Initial 
rotation of massive stars is not yet fully constrained, angu¬ 
lar momentum gain during accretion in massive binaries is 
not fully understood and angular momentum loss with stel¬ 
lar winds for massive stars is still uncertain (see Sec. |4] for a 
more detailed discussion). Therefore, we study a broad spec¬ 
trum of massive stars with BH, to which we refer as B-BH 
binaries and which are defined below. Some of these binaries 
may produce or are born with Be or Oe star. 

In our simulation B-BH systems are binaries satisfying 
a following criterion: 

BH + MS binary 

10 < Porb < 1200d (3) 

3 < MS mass < 30 Mq. 

Note that this specific choice of orbital periods along with 
MS secondaries leads to population of wind-fed binaries. We 
have chosen the orbital periods such, that they are within 
the expected range for Be X-ray binaries (see Sec. 11.51) . 
Also note that MS mass range includes not only B stars 
(< 15Mq) but also some O stars as to corres pond to the ob¬ 
served Be/Oe population dNeeueruelal (Il998ll and references 
therein). 

We also define a subpopulation of B-BH binaries 
with properti e s sim ilar to these observed for MWC 656 
ICasares et alJ ll2014l ~l as: 

BH + MS binary 
55 < Porb < 65d 

3.8 < Mbh < (^.QMq (4) 

10 < Mms < I6M0 
e < 0.5 

The masses h ave been chose n with in the one sigma limits 
determined bv ICasares et al] ll2014ll : the eccentricity range 
was established in a private communication with the team 
which discovered the first BH-|-Be binary; the range for the 
orbital period was chosen arbitrarily. 

Note, that to obtain the information on Be-BH systems 
or systems similar to MWC 656 containing a Be star, all 
numbers/rates we get for the groups defined above need to 
be multiplied by a factor fse, which is not well known. If 
the fraction of systems where the B star shows the Be phe¬ 
nomenon were similar to the observed fraction in Milky Way 
and Magellanic Clouds clusters, the reduction factor /se 
would be about 0.3 (see Sec. 


3 RESULTS 

In the following text we refer to initially more massive star 
as primary (BH progenitor), and initially less massive star 
(MS star that is potential Be or Oe object) as secondary. 

3.1 Overall properties of B-BH binaries 

The total number of B-BH systems formed in our standard 
model simulations over entire 10 Gyr of evolution of the 



Figure 1. Orbital period distribution for overall group of B-BH 
binaries (red solid line) and MWC 656-like subpopulation (blue 
dashed line). Predicted current Galactic populations are shown. 
Note that overall number of B-BH binaries (39) is much larger 
than for MWC 656 - like systems (0.01 - it is a probability of 
having one MWC 656 - like system). 


Galactic disk is = 8, 700, while number of MWG 656- 

like systems is A^Mwcese = 13. 

Note that our predictions for MWG 656-like systems are 
subject to errors from small number statistics. But these 
errors are smaller than ones associated with evolutionary 
uncertainties. Various systems that we generate have a range 
of lifetimes during B-BH (with an average of 45 Myr) 
or MWC 656 (with an average of ~ 6 Myr) stage. We use 
the lifetime of each system to assess the probability that it 
is present in the current Galactic disk population. We use 
N = 10'^ of Monte Carlo realizations of the formation time 
(drawn from uniform distribution) of a given system in range 
0 — 10 Gyr and check how many times n a given system is 
present at the current Galactic age (10 Gyr). The current 
number of given binary population is given by n/N. This 
number may be smaller than one and then it indicates how 
low is the probability of one system to exist in the current 
predicted population. 

The total number of B-BH systems found in our simu¬ 
lations to be present at the current moment in the Galactic 
disk is = 39, while number of MWG 656-like systems 

is well below one: = 0.007 (i.e., the probability of 

having one system at present is ~ 1%). 

In Figures [T][^ [3] and [4] we show distributions of orbital 
period, black hole mass, companion mass and orbital eccen¬ 
tricity for B-BH binaries (defined by Eq. [3]). We also show a 
subpopulation of B-BH binaries that resemble the observed 
properties of MWG 656 (defined by Eq. [3|. In Figures [S] 
and [6] we show the two-dimensional distributions of orbital 
period vs. B star mass and eccentricity vs. BH mass for B- 
BH binaries. With a white rectangle we mark the region 
corresponding to systems defined as MWC 656-like (Eq. |4]). 

In Fig. [7| we present the sensitivity of total number of 
MWC 656-like systems formed over entire 10 Gyr in the 
Galaxy on the adopted orbital period range in the systems’ 
definition Eq. [d] 

There is a two-stage evolution leading to the formation 
of B-BH binaries; GE followed by core-collapse/SN (forms 
BH). In the first stage the massive primary expands after 
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Figure 2. Black hole mass distribution for overall group of B-BH 
binaries (red solid line) and MWC 656 - like subpopulation (blue 
dashed line). Predicted current Galactic populations are shown. 



M2 [MJ 

Figure 3. B star mass distribution for overall group of B-BH 
binaries (red solid line) and MWC 656 - like subpopulation (blue 
dashed line). Predicted current Galactic populations are shown. 



Figure 4. Eccentricity distribution for overall group of B-BH 
binaries (red solid line) and MWC 656 - like subpopulation (blue 
dashed line). Predicted current Galactic populations are shown. 
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Figure 5. Two-dimensional distribution of orbital period (x 
axis) and B star mass {y axis) for B-BH binaries. The numbers 
next to the color bar indicate the logarithm of the predicted cur¬ 
rent Galactic B-BH population. The white rectangular box indi¬ 
cates the area of MWC 656-like systems defined in Eq. [4| 
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e 

Figure 6. Two-dimensional distribution of eccentricity (x axis) 
and BH mass (y axis) for B-BH binaries. The numbers next to the 
color bar indicate the logarithm of the predicted current Galactic 
B-BH population. The white rectangular box indicates the area 
of MWC 656-like systems defined in Eq. [4| 


MS and initiates RLOF. Due to the high mass ratio (typ¬ 
ically Mprimary/Msecondary ^ 3) mass transfer is unstable 
and leads to CE. At the onset of CE, our massive primary 
(^ 30-60M©) is most likely (~ 90%) in the Core Helium 
burning (CHeB) phase and has a well developed convective 
envelope. Only in a relatively small number of cases (~ 10%) 
the primary is on the Hertzsprung gap (HG) with either a 
radiative or shallow convective envelope^. In the case of 


2 


Note that these numbers have changed significantly since 
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Figure 7. Sensitivity of formation efficiency of MWC 656-like 
systems (as defined in Eq. 0. For our study we have adopted 
definition 55 < Porb < 65d for MWC 656-like systems in respect 
to orbital period. For this definition we produce 13 systems in 
entire 10 Gyr of Galaxy evolution (see the blue dot). If we widen 
the period range by a given number of days (horizontal coordi¬ 
nate) on both sides of the measured orbital period of MWC 656 
= 60.37 d) then we obtain specific number of MWC 656-like 
systems (vertical coordinate). 


HG donor it is not at all clear whether CE develops at all 
even for high mass ratios (unpublished MESA simulations). 
This 90-10% division is an evolutionary selection effect. The 
secondary star is typically of much lower mass (--^ 5-15Mq; 
see Fig. El than primary. Therefore, it is easier for a binary 
to survive the CE phase for small CE envelope mass. HG 
stars have massive and tightly bound envelopes, while the 
envelope mass and binding energy decreases rapidly during 
CHeB phase. Both as a result of expansion and for high 
metallicity also due to intense wind mass loss. After CE 
initially wide systems evolve to much shorter orbital peri¬ 
ods 10-100 day; see Fig. [!}. Massive primaries lose their 
entire H-rich envelope, and become Wolf-Rayet stars. We 
assume that CE is very rapid and that the MS secondary 
does not have time to accrete. In our simulations we allow 
only compact objects (NS and BH) to accrete a small but 
significant amount of mass during CE. Note that at the mo- 
men t the CE phase is fa r from being understood (see how¬ 
ever [iranoTO^eT^y [201^ and our assumptions on accretion 
physics in CE are subject to verification. 


iBelczvnski &: Ziolkowsldl ll2009ll . who reported only ~ 4% of pro¬ 
genitors of B-BH binaries to go through CE with a donor on 
CHeB stage. This apparent discrepancy is a direct consequence 
of updates of input physics made in StarTrack code. It mostly 
originates from change of common envelope physics. In the up¬ 
dated code the iA parameter that describes the binding energy 
of the donor star is calculated based on stellar radius, metallicity 
and evolutionary stage. In the previous study the authors used a 
constant value A = 1, which is now replaced with a typical phys- 
ical value of ~ 0.1 — 0.2 for BH progenitors on HG (see Fig. 3 in 
iDominik et al.l ll20l2tl 1. It means, that now A is too small to allow 
HG donor to survive CE, and these systems now are found to 
typically merge during CE. The smaller lambda values preferen¬ 
tially select donors with small envelope mass to survive CE. For 
massive BH progenitors stars lose their H-rich envelopes during 
CHeB and that is what we note in our new results. 
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Figure 8. Typical evolution that may lead to the formation of 
MWC656 - like system. For the detailed description of the evo¬ 
lutionary history see Sec. I3.2I Note that the two most important 
evolutionary factors are common envelope phase that brings or¬ 
bital period close to the current observed value and supernova 
Ib/c that forms rather light black hole similar to the one residing 
in MWC656. 

The notation: ZAMS - zero age main sequence, MS - main se¬ 
quence, HG - Hertzsprung gap, CHeB - core helium burning, 
HeMS - helium main sequence, BH - black hole, CE - common 
envelope, SN - super nova 


In the second stage, the Wolf-Rayet star explodes in 
Type Ib/c SN and forms a BH. Alternatively, for very mas¬ 
sive stars we assume direct BH formation without accompa¬ 
nying SN. Our models for single stars allow for Wolf-Rayet 
star formation above Mzams ~ 20Mq and for direct BH 
above Mzams ~ 4OM0. Roche lobe overflow episodes in bi¬ 
naries may signihcantly shift these boundaries up or down 
depending on a given binary conhguration and evolutionary 
stage of its components. The associated with SN nentrino 
losses and potential mass loss and natal kick affect the or¬ 
bit (altering eccentricity and the semi-major axis). A BH 
binary with MS companion is formed. Such a binary may 
either fall right within our criteria for B-BH or MWC 656- 
like object directly, or it may evolve to satisfy these criteria 
at a later stage. The only process altering the binary orbit 
at this stage is wind mass loss from the MS star increasing 
the orbital separation (but this effect is very small for the 
MS stars in the mass range we consider). For close binaries 
tidal interactions (synchronizing MS stars and circularizing 
the orbit) also play a role. Since we have chosen our lower 
limit on the orbital period to be rather large (10 days), the 
systems that are subject to efficient tidal interactions are 
not typical progenitors of B-BH/MWC 656-like objects. 


3.2 The formation of MWC 656 - like system 

Here, we present a typical evolutionary scenario that leads 
to the formation of a MWC 656 - like system (see Fig. EJ. 

We start the binary evolution with two components on 
the Zero Age Main Sequence (ZAMS) with Mi = dlM© 
(primary) and M 2 = IAMq (secondary). The initial semi¬ 
major axis of the orbit is a = 5.3-10®R© and its eccentricity 
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is e = 0.42 (the orbital period is Porb = 6000 day). After 
4.8 Myr the primary with a mass Mi = 35Mq hnishes core 
hydrogen-burning and enters the HG. At this phase the pri¬ 
mary significantly expands and tidal forces begin to circu¬ 
larize the orbit. Next, the primary enters CHeB expanding 
toward its Roche lobe. The orbit becomes fully circularized. 
This phase ends with the primary (Mi = 17M©) overfilling 
its Roche lobe and initiating the common envelope phase. 
At t = 5.1 Myr a close binary emerges out of the CE - a 
Ml = ISM© helium core of the primary on a relatively close 
orbit (separation decreases from a = 5300i?© to a = 1427?© 
corresponding to decrease in orbital period from Porb = 8000 
day to Porb = 38 day) around the mostly unaffected main se¬ 
quence secondary M 2 = 14M©. After ~ 0.35 Myr the helium 
star primary finishes nuclear burning and its mass decreases 
to Ml = lOM© due to strong Wolf-Rayet type winds. Just 
before the supernova (SN) explosion the semi-major axis of 
the orbit is a = 167P© and the orbital period Porb = 52 
day (expansion due to the wind mass loss). The primary 
explodes in a type Ib/c SN (ejected mass ~ 2.5M©) and 
forms a light black hole with mass Mi = 5.6M© (we have 
assumed 10% mass loss in neutrino emission). We obtain 
a natal kick from an asymmetric mass ejection scenario and 
for this particular system the magnitude of the kick is 130 
kms“^ in such direction that it changes the separation to 
a = 165P© (Porb = 56 day) and eccentricity to e = 0.01. 
Depending on the orientation of the kick, the post-SN ec¬ 
centricity may vary in a broad range. For small post-SN 
eccentricities (e < 0 . 1 ) systems have a likely chance to form 
BH-NS binary at the end of evolution, while more eccentric 
systems tend to finish evolution in CE mergers (see Sec. 13.51) . 
Over the next 9 Myr this relatively close binary remains 
almost unchanged with slight orbital expansion due to wind 
mass loss from the main sequence secondary (a = I 68 P 0 , 
Porb = 58 day). Throughout this phase the system is a wind- 
fed X-ray binary and it meets our criteria Eq. 0] and we tag 
it as the MWC 656-like system. 

Initially MWC 656 progenitors are two massive stars. If 
placed on a short orbit, such stars would begin interacting 
while on the MS or at the beginning of the HG and such a 
RLOF would mo st likely lead either to component merger 
llSana et al.ll201 or would deplete the primary mass below 
the threshold of BH formation. For wide binaries, it takes 
a CE phase to decrease the orbital size to the currently ob¬ 
served period of MWC 656 (as described in our example). 
For very wide systems either CE is never encountered or if it 
is there is only very little mass in CHeB donor envelope and 
orbital decrease is not efficient enough to produce orbital 
period observed for MWC 656. 

3.3 Parameter study 

Our results are subject to a number of evolutionary uncer¬ 
tainties. We have performed several additional calculations 
that probe the most obvious uncertainties in the formation 
of B-BH binaries. In models Vi and V 2 we have altered the 
energy balance in CE evolution that is required in all chan¬ 
nels leading to B-BH binary formation. In models V 3 and 14 
we have changed our treatment of BH formation (natal kicks 
and BH mass). These changes lead to different formation ef¬ 
ficiency of B-BH binaries. In Table [T] we list total number of 
B-BH binaries formed over entire 10 Gyr of Galactic evolu¬ 


tion as well as their current predicted number in Galaxy for 
all models. The same numbers are listed for MWC 656-like 
systems. 

The current number of B-BH systems is predicted at 
the level of ~ 10-100. This number is rather sensitive to the 
adopted assumption on CE efficiency. For model in which 
we allow only 10 % {acE = 0 . 1 ; model 14) of orbital en¬ 
ergy to be used in CE ejection we predict only 7 B-BH 
systems to be currently present in our Galaxy. For much 
more efficient ejection, with 5 times of orbital energy is 
used {acE = 5.0; model V 2 ) for CE ejection we find many 
more B-BH systems: 131. We only allow orbital energy to 
increase in such arbitrary way as to mimic the possibility 
that bounding energy of CE is much lower than used in 
current predictions. We use bindi ng energy estimates from 
IXu fc Lil ll2010bl ') dXu k, Lill2010al ~l that give binding energy 
of envelope for a given star radius, metallicity and star evo¬ 
lutionary stage. However, the envelope internal energy may 
lead to much easier ejection and it was estimated that it 
is pote ntially realistic to decrease bi nding energy by factor 
of 5 lllvanova fc Chaichenet j|201ll ). That is what we have 
employed in model V 2 . For our reference model we have used 
100% of orbital energy for CE ejection [acE = 1.0; model 
So). 

Typical B-BH formation starts with a rather massive 
star (BH progenitor) that forms massive envelope after main 
sequence. Ejection of massive CE (10 — OOM©; at the end 
of Hertzsprung gap) by a typical B star (~ 5-15M©; see 
Fig. O is rather hard. In particular, for low ejection effi¬ 
ciency (acE = 0.1) it leads to a decrease in number B-BH 
binaries. For model in which we allow for lowered binding 
energy (or high acE = 5.0) the B-BH number increases. Al¬ 
ternatively, binary channels are naturally selected in which 
the envelope of a massive primary is depleted by evolution 
(winds and core growth) and then late case C RLOF leads to 
CE development (with CHeB donor) and B-BH formation. 

The delay ed SN mode l allow s for the formation of low 
mass BHs isee lFrver et al.l (l2012l H. as opposed to our stan¬ 
dard model in which we do not allow formation of BHs in 
the mass gap: 3 — 5M© llBelczvnski et alll2012ll . We note 
no significant change in number of B-BH binaries between 
these two models (standard vs. delayed SN). There is a small 
increase of B-BH binaries in the model with no BH kicks, 
since in this case some progenitors are not disrupted upon 
BH formation. Since the majority of BHs in B-BH binaries 
are predicted to have rather high mass (peak of mass dis¬ 
tribution at 7 — 8 M©; see Fig. n then they receive small or 
no natal kicks in our standard model. Therefore, the change 
to zero BH kick has a small effect on the overall B-BH pop¬ 
ulation. This is quite different for MWC 656-like systems, 
for which we note factor of ~ 10 increase of current Galac¬ 
tic number in model V3. As low-mass BHs (as observed in 
MWC 656) receive non-zero kicks (binary disruptions) at 
the formation in our standard model, thus in model V3 (no 
kicks) we note significant increase of MWC 656-like systems 
in the overall population of B-BH binaries. 

We note that MWC 656 eccentricity and peculiar space 
velocity are consistent with no or small natal kick. The high 
BH natal kicks (above 100-150 kms“^) may be excluded by 
the analysis presented in the Appendix. In our particular ex¬ 
ample of the formation of MWC 656-like systems (presented 
in Fig[ 8 ]) the BH is formed with a moderate 3D natal kick 
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Table 1. B star + black hole binary numbers in Milky Way“ 


model 

MWC 656-like 

B-BH 

comment 

So 

13 

(0.007) 

8,700 

(39) 

standard evolution 

Vi 

0 

(0) 

1,600 

(7) 

"CE = 0.1 

Vi 

34 

(0.026) 

55,800 

(131) 

«CB = 5 

Vs 

71 

(0.063) 

12,800 

(63) 

no BH kicks 

Vi 

8 

(0.004) 

13,300 

(50) 

delayed SN engine 


“ The total number of B-BH and MWC 656 - like systems formed in simulations 
over 10 Gyr of evolution of the Galaxy; in parenthesis we list the current number of 
systems in the Galaxy. 


Vkick ~ 130 kms“^ and mass loss of ^ SM©. This natal kick 
was oriented in such a way that eccentricity was essentially 
not affected (cpostSN = 0 . 01 ), but the systemic velocity was 
increased by 50 kms“^. This is close to the values presented 
in the Appendix: e < 0.14 and KpaceMwcess < 37.4 kms“^. 
Within our population of MWC 656-like systems we find bi¬ 
naries that are fully consistent (in a 1 -sigma range) with the 
Appendix estimates: (i) for natal kick of 23 kms“^ and mass 
loss of ~ 3 M 0 we get e = 0.1 and KpaceMwcese = 22 kms“^ 
and (ii) for natal kick of 80 kms“^ and mass loss of ~ ZMq 
we get e = 0.12 and l^paceMwcsss = 36 kms“^. The above 
examples were obtained with our standard model, which as¬ 
sumes non-zero BH (but smaller than NS) kicks. Our model 
for non-zero kicks is based on the asymmetric mass ejection 
from pre-supernova star. The amount of mass loss (that sets 
the natal kick value) i s base d on the supernova models pre¬ 
sented in iFrver et al.l (120121 '). Apparently these models can 
explain both the eccentricity and space velocity of MWC 
656. Obviously, our models with no BH natal kicks are also 
fully consistent with the eccentricity and peculiar space ve¬ 
locity estimates. Therefore, we cannot distinguish between 
the two models (asymmetric mass ejection versus no natal 
BH kicks). 


3.4 On the origin of the Be phenomenon in 
Be-BH systems 

The emission lines of regular main sequence Be stars are gen¬ 
erally considered to result from an outflowing disk of a star 
rotating at or at a substantial f raction of its Kepler ian ro¬ 
tation rate or break-up velocity llRivinius et al.ll2013lj . This 
raises the question about the origin of the spin of the Be 
star. This may either be reflecting the birth spin of the star 
or it may be the consequence of interaction in the binary 
system. 

The first case is plausible, since young massive 
stars are obser ved to rotate at a wide range of rota¬ 
tion rates (e.g. lAbt et al.l l2002l : iRamfrez-Agudelol l2013l : 
ISimon-Dfaz fc Herrerd 2014l l. It does however require the 
Be star to retain its spin through all phases of binary inter¬ 
action including the common envelope. This phase is poorly 
understood. Detailed stellar structure models show that ro¬ 
tating main sequence stars naturally tend to spin up their 
outer layers towards the Keplerian rotation rate. This does 
however require that the outer layers are well coupled with 
the contracting core. It also requires that angular momen¬ 
tum loss in the form of stellar winds is not significant, which 
is the case for most B-type main sequence stars and at low 


metallicity for the later O t ype stars as well (lEkstrom et al.l 
l2008bl : Ide Mink et al.l[2013lj . In the currently favored picture 
this can cause stars that are rotating sufficiently fast to reach 
break-up as they evolve and show the emission phenomenon 
until they leave the main sequence. 

The second case requires a phase of mass transfer from 
the progenitor o f the current black hole to the B star. 

I Pols et ffi] lll99lll showed that spin up by mass transfer is 
very efficient: accreting just a few percent of its mass from 
an accretion di sk is suffic i ent to bri ng a star to break-up . 
Simulations bv iPols et al.l lll99ll j and Ide Mink et ^ (l2013lj 
show that, given the high binary frequency among massive 
stars, this is expected to be a very important or possibly 
even dominan t channel for the formation of rapidly rotating 
stars (see also lShao fc Li|[2014h . 

In our models we find that none of our MWC 656-like 
systems (containing both B and Be stars) formed through 
a formation channel that included a phase of stable Roche 
Lobe overflow. So our simulations give preference to the first 
case, where the spin of the B star is the spin resulting from 
birth. However, it must be noted that in our model we rely 
on simplified prescriptions of the stellar structure and there¬ 
fore an approximated treatment of the response of a star 
to Roche lobe overflow. One possibility is spin up during a 
(presumably short phase of) stable mass transfer preceding 
the CE. Another interesting alternative is stable mass trans- 
fer through atmo spheric or wind Roche lobe overflow (e.g. 
lAbate et al.ll 20 f^ . 


3.5 The future evolution of MWC 656-like system 

The fate of MWC 656-like system strongly depends on the 
secondary mass on ZAMS. We selected two broad evolu¬ 
tionary categories; binaries with relatively low-mass 10 < 
M2zams < 13M0 and high-mass 13 < M2zams < I 6 M 0 
secondaries. The summary of future evolution of MWC 656- 
like systems and the estimates of BH-NS formation chances 
are given in Table [2] 


3.5.1 Progenitors of wide BH-NS binaries 

Here we describe, the future evolution of MWC 656-like sys¬ 
tems with the initial (ZAMS) secondary mass smaller then 
13M0. This group consists of ~ 38% of all MWC 656-like 
systems formed in Galactic disk, with 15.9% forming wide 
BH-NS systems and 22.6% are being disrupted in second 
SN. 

We start with the binary consisting of a BH with mass 
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Table 2. Future evolution of MWC 656-like binaries® 


Channel 

/form 

Evolutionary history^ 

Mergers^ 

CE/RLOF 

Fate‘s (BH-NS): 

Close Wide Disrupted 

B-BH:la 

15.4% 

CE1(4-1) SNl MT2(14-2) MT2(14-9) ECSN2 

0 % 

0 % 

15.4% 

0% 

B-BH:lb 

23.1% 

CE1(4-1) SNl MT2(14-2) SN2 

0 % 

0 % 

0.5% 

22.6% 

B-BH:2a 

7.7% 

CE1(4-1) SNl CE2(14-4) MT2(14-7) SN2 

0 % 

5 . 6 % 

0.7% 

1.4% 

B-BH:2b 

53.8% 

CE1(4-1) SNl CE2(14-2) MT2(14-7) SN2 

38.4% 

10 . 7 % 

1.3% 

3.4% 


® We list only formation channels of MWC 656-like systems which are defined by Eq. 

^ Sequences of different evolutionary stages: CEl and CE2: common envelope with a primary and secondary as a 
donor, respectively; MT2: non-conservative mass transfer with a secondary as a donor; SNl and SN2: type Ib/c 
supernova of the primary (black hole formation) and secondary (neutron star formation), respectively; ECSN2: 
electron capture supernova of secondary (neutron star formation). 

Numbers in parenthesis denote evolutionary stage of primary—secondary: 1 - main sequence, 2 - Hertzsprung gap, 4 
- core helium burning, 7 - helium main sequence, 9 - helium giant branch, 13 - neutron star, 14 - black hole. 

® This is probability that two binary components merge in RLOF or CE events that are encountered between the 
two SNe events. 

Outcome of future evolution of MWC 656-like systems; close (delay time from ZAMS to BH and NS merger 
shorter than 10 Gyr) or wide BH—NS systems or disrupted BH and NS objects may form. 


Ml = 5.7Mq, and a secondary star that has just entered 
HG with mass M2 = IO.2M0 (so right after MWC 656-like 
phase) and with separation a = 1537?© (Porb = 52.2 day) 
and eccentricity e = 0.28. 

The secondary quickly expands while crossing HG and 
it initiates stable Roche lobe overflow (RLOF) onto the 
BH. Rapid expansion does not allow for tidal circulariza¬ 
tion. We circularize the system instantaneously at the on¬ 
set of RLOF, we take the periastron distance as the new 
separation (a = HOP©, Porb = 34 day) of circular or¬ 
bit (e = 0). The accretion onto BH is Eddington limited, 
and the rest of mass leaves the system with the BH spe¬ 
cific angular momentum. During RLOF the separation first 
decreases to a = lOlP© (Porb = 32 day) and then in¬ 
creases to a = 364P0 (Porb = 263 day). The binary com¬ 
ponents go through a mass ratio reversal, with BH mass 
Ml = 7.IM0 and the secondary mass M2 = 2.3Mq at the 
end of RLOF. The secondary becomes a naked helium star 
with core Helium-burning that lasts about 4 Myr. The low 
mass helium secondary begins to significantly expand after 
it becomes an evolved helium star (Helium shell burning) 
and it initiates another episode of RLOF. At this point, the 
secondary mass decreases to M2 = 2.IM0 and the binary 
separations increases to a = 423P0 (Porb = 333 day). As 
a result of the RLOF the BH mass reaches Mi = 7.2Mq 
and secondary mass is depleted to M2 = I.6M0, while the 
orbit expands to a = 623P0 (Porb = 608 day). The sec¬ 
ondary ends its evolution as a neutron star (NS) of a mass 
M2 = I.26M0 created in electron capture SN. All binary 
systems survive the explosion as we assume no natal kick 
for electron capture SN and very little mass was lost in the 
process. This evolutionary channel is marked as “B-BH:la” 
in Table [2] We note the formation of wide BH-NS system 
with chirp mass of Mo = 2.4M0 and a very long merger 
time tmerger = 3.9 X 10® Gyr. 

For slightly more massive secondaries iM2zAMs ~ 
12 Mq) than in the above case, the second RLOF episode 
is avoided and a NS with mass M2 = I.IM© is created in 
a type Ib/c SN explosion. The chance of survival of such a 
binary is only ~ 2% due to frequent natal kick disruptions. 
The systems that survive SN explosion form wide BH-NS 


binaries with merger times exceeding the Hubble time. This 
evolutionary channel is marked as “B-BH:lb” in Table [21 

3.5.2 Progenitors of close BH-NS binaries 

In this section we describe the fate of MWC 656-like system 
with a secondary mass on ZAMS larger than 13Mq. Such 
systems make up 62% of the population of MWC 656- 
like binaries formed in Galactic disk. Although a significant 
fraction of these systems merge in ensuing CE (38.4%), and 
some fraction gets either disrupted in the second supernova 
explosion (4.8%) or form wide BH-NS systems (2.0%), a 
sizable fraction forms close BH-NS systems (16.3%). 

We start with the binary consisting of a BH with mass 
Ml = 5.7M0, and a secondary star that just has entered 
HG with mass M2 = 13.5M0 (right after MWC 656-like 
phase) and with separation a = I68P0 (Porb = 58d) and 
eccentricity e = 0.01. 

Typically, lower mass secondaries {M2zams ^ 14.5M0) 
evolve through HG and enter CHeB and then initiate CE 
(formation channel “B-BH:2a”; see Table[2]). The chances of 
survival of this CE are close to unity due to small CE en¬ 
velope mass. Higher mass secondaries {M2zams ^ 14.5M0) 
initiate CE while still on HG (“B-BH:2b”). The chances of 
survival of this CE phase are only about one third due to 
large CE mass. It is possible that instead of CE some of these 
systems evolve through fast (on thermal timescale) but sta¬ 
ble RLOF. In such a case these channel will contribute very 
little (if any) to the formation of close BH-NS binaries. 

After CE, the secondary loses most of its mass and be¬ 
comes a naked helium star with mass M2 = 3.3Mq. The 
orbit is circularized and decreases in size to a = I.5P0 
(Porb = L7 hour). A very compact binary is formed. The 
core-Helium burning secondary expands and initiates a sta¬ 
ble RLOF. After 2 Myr the secondary becomes an evolved 
helium star. The RLOF continues. The BH mass increases 
to Ml = 6M0, while the secondary mass decreases to 
M2 = I.6M0 and at the time of the SN explosion the or¬ 
bit has expanded to a = dP© (Porb = 8.5 hour). After 16.5 
Myr from the beginning of the evolution at ZAMS, the sec¬ 
ondary explodes in SN type Ib/c, forming a NS of a mass 
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M 2 = I.IMq. Chances for a binary survival are very high 
as these systems are very compact. 

For one particular natal kick, the post-SN orbit becomes 
eccentric e = 0.46 and the semi-major axis increases to 
a — 5.7Rq (Porb = 14.4 hour). The close BH-NS binary 
is formed with the chirp mass Me = 2.IM0 and merger 
time fmerger = 1.5 Gyr. 


3.6 Future evolution of MWC 656 


Based on the results presented in the above sections we 
can describe the future evolution of a binary that resem¬ 
bles MWC 656 (see Fig. [5]). We start the evolution of the 
binary with a black hole mass Mi = 5.35M©, B star mass 
M2 = 13.5M©, semi-major axis a = 172P© and eccentricity 
e = 0.1 and this resu lts in an orbital period of Porb = 60.37d 
dCasares et al ][2oT3). The parameters are all within the ob¬ 
servational errors of the original observational estimates and 
allow for the system to evolve along our evolutionary chan¬ 
nel “B-BH:2a”. The massive primary (Mzams = 30 — 35Mq) 
took ~ 6.2 Myr to form a BH. A secondary MS lifetime 
is ~ 14.8 Myr. This gives us a typical lifetime of the B- 
BH binary phase of tMWC656 = 8.6 Myr. After finishing MS 
evolution, the secondary enters HG and starts to burn Fle- 
lium in its core. During this stage a CE phase is initiated. 
The orbit significantly decreases after CE and the secondary 
becomes a low-mass naked helium star. The secondary ex¬ 
pands again while core-helium burning and initiates stable 
Roche lobe overflow. A BH increases its mass during CE 
(by ~ O.I2M0) and during stable RLOF (by ~ 0.3Mq) to 
Ml = 5.8M0. Finally, after tevoi = 17.2 Myr from ZAMS 
the secondary explodes in type Ib/c supernova and forms 
a NS with mass M2 = I.IM©. SN neutrino and mass loss 
(~ O.32M0) and natal kick (drawn from Maxwellian with 
a = 265 km s“^) may lead to a disruption of a binary (prob¬ 
ability of 13.2%), formation of wide non-coalescing BH-NS 
system (10.1%) and formation of close: coalescing within 10 
Gyr BH-NS system (76.7%). There is very high probability 
of close BH-NS formation /dose = 0.77 as the circular pre- 
SN orbit is very compact and hard to disrupt (a = 4.OP©; 
Porb = 8.6h). We use a high number of Monte Carlo experi¬ 
ments to assess these probabilities and we show the resulting 
distribution of delay times in Fig. HOI The close BH-NS sys¬ 
tems that may form out of MWG 656 have average (and 
median) delay time of 1.9 Gyr (and 0.9 Gyr). 

Note that initial stages of our prediction for future evo¬ 
lution of MWG 656 (Fig[5)) resem ble a recently disc overed 
ULX source: P13 in NGG 7793 (iMotch et al. [20l3). P13 
is a binary with a 3 — IbM© BH and rsj ig — 23Mq B9Ia 
companion. It was estimated that the system is experienc¬ 
ing the supercritical RLOF. The companion is proposed to 
be at the end of MS or at the beginning of HG. The or¬ 
bital period is about 64 days. The only significant difference 
between P13 and MWC 656 is the mass of BH companion 
star. Since MWC 656 star has mass smaller than that in 
PG13, the RLOF will start later in its evolution. It is pre¬ 
dicted that the onset of RLOF will begin when the star in 
MWC 656 will already finish crossing HG and will start Gore 
He burning. CE may be preceded by the short and stable 
high mass transfer rate phase (then the MWC 656 would re¬ 
semble P13). However, it seems that due to the existence of 


time [Myr] 


6.2 

5.35Mg 

BH 

• • 

13.5Mg 

MS 

a = i72i2g 

P=60d 

e= 0.1 

14.7 

5.35Mg 

BH 


13.2Mo 

HG 

a = 175i?0 

P=62d 

e= 0.1 

14.8 

5.35Mq 

BH 

m 

i 13.2M^ 

i CHeB 
^ before CE 

0 = 175 
P=e3d 
e= 0.1 

14.8 

5.47M^ 

BH 

•• 

3.2M^ 
HeMS 
after CE 

a=1.4Ro 

P=1.6/i 

e= 0.0 

15.2 

5.47Mg 

BH 


3.2M^ 

HeMS 
before MT 

a=lAR^ 

P=1.6h 

e=0.0 

16.1 

5.52Mg 

BH 

•• 

2.9Ma 
HeMS 
after MT 

a=lSR^ 

P=1.8h 

e= 0.0 

17.1 

BH 


2 . 8 M 3 
evolved He 
before MT 

a= 1 . 6 i?o 

P=2h 

e= 0.0 

17.2 

5 . 8 M 3 

BH 

•• 

l.SMg 
evolved He 
after MT 

c,=4.0i?^ 

P= 8 . 6 /i 

e= 0.0 

17.2 

5.8M^ 

BH 

• • 

I.IM 3 

SN - NS 

o = 14R 
P=2.2d 
e= 0.8 


Figure 9. The future evolution of MWC 656 system (see 
Sec. l3.6ll . The binary evolves through CE and stable RLOF phase. 
If the binary survives CE phase the secondary star will form a 
light NS in supernova explosion (f.lM©). Since CE significantly 
decreases the orbital separation, the binary is very likely to sur¬ 
vive supernova mass loss and natal kick and to form a close BH- 
NS system (probability of 77%). 

We use the same notation as in Fig[Sl with the addition of MT - 
mass transfer, NS - neutron star. 


deep convective envelope that forms during CHeB and quite 
high mass ratio (2.5) the development of CE is very likely. 

Another object, SS433, was pro posed to be a mas sive 
star engulfing a BH in its envelope llClark et al.ll2007h . It 
may be the only known case of a massive binary undergo¬ 
ing CE. So far all the other bina ry mergers/CE events are 
restricted to low-mass stars fe.g.. Il^chanek et al.l ll2014ll L 


3.7 Empirical LIGO/VIRGO detection rates 

The lifetime of MWC 656 may be estimated from our evo¬ 
lutionary calculations presented in Sec. 13.61 Imwcsss = 8.6 
Myr. If we assume that only one such system is present cur¬ 
rently in the Galaxy we obtain the Galactic birth rate of 
7?.birth ~ l/tMWC656 Under assumptions that star formation 
in Galaxy was constant and nothing special or extraordi¬ 
nary was required to form MWC 656. We have shown in 
Sec. 13.61 that the probability of forming a close BH-NS bi¬ 
nary out of MWC 656 is /©ose = 0.77. Since the delay times 
are relatively short (median of the distribution is 0.9 Gyr; 
see Fig. HOJ as compared with the Galactic disk age (10 Gyr) 
we can estimate the Galactic merger rate as 

7^MW = /close7^bi^th = , = 0.089 Myr’^ (5) 

iMwcese 

We have converted the Galactic merger rates to the ad¬ 
vanced LIGO/Virgo detection rates (TT-ligo) assuming the 
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Figure 10. The delay time distribution for BH-NS systems 
formed out of MWC 656-like population defined by Eq. U] (red 
solid curve; average 1.8 Gyr and median 0.8 Gyr) and for BH-NS 
systems formed from the system with exact observed MWG 656 
properties as defined in Sec. 13.61 (blue dashed curve; average 2.0 
Gyr and median 1.1 Gyr). 


constant density of Milky Way-like galaxies at the level 
Pgai = 0.01 Mpc~® in local Universe. We have adopted 
do = 450 Mpc as the advanced LIGO/VIRGO horizon 
for NS-NS binary (optimally oriented source with signal-to- 
noise ratio of 8) with chirp mass Ale,nans = {Mi + 

= 1.2M0 where individual NS masses are Mi = 
M 2 = IAMq. The horizon for a double compact object 
with a given chirp mass A4c,dco is calculated with d = 
do(AIc,dco/A4c,nsns)®^®. Finally, the detection rate is ob¬ 
tained with: 


Tcligo = Pgai — 



k Alcjdco 

\ AIc,nsns 


15/6 

TImw 


( 6 ) 


where factor /pos = 2.26 takes into account the non-uniform 
pattern of detector sensitivity and random sky orientation of 
sources. For our case of BH with Mi = 5.8Mq and NS with 
mass M2 = I.IM0 we have Mfc.dco = 2.IM0 and the cor¬ 
responding detection rate of 7?.ligo = 0.115 yr“^. In other 
words, the existence of MWC 656 binary implies 1 BH-NS 
advanced LIGO/Virgo detection in 9 years. 

This estimate is subject to a number of uncertainties. 
For example, if we lower mass of the B star ~ IOMq then 
such a system will form only wide BH-NS binaries (see 
evolutionary channel “B-BH:la” in Table II) and we get 
'Ti-iAGO = 0 yr“^. If we lower B star mass to 11 — I2M0 
such system will most likely get disrupted by core collapse 
SN (see evolutionary channel “B-BH:lb” in Table[2)) and we 
also obtain 7?.ligo = 0 yr“^. These null rates are the di¬ 
rect result of our criterion on development of CE. For lower 
mass stars, instead of CE we encounter stable RLOF and 
therefore we form much wider systems that are either dis¬ 
rupted or do not merge within Hubble time. It is worth 
noting that the CE development and its inner mecha nism 
is still far from being understood lllvanova et al. I l2013l) . On 
the other hand we may increase the rate by increasing the 
B star mass to IhM© (this shortens the B lifetime to 
tMWG656 = 6.4 Myr) and lowers the kick velocities by factor 
of 2 (this increases chance of SN survival to /dose = 0.92) to 
obtain 7?.ligo = 0.187 yr“^ (1 detection in 5 years). Lower 


NS kicks in interacting binary systems are supported for ex¬ 
ample b y the observed ratio of single to binary millisecond 
pulsars llBelczvnski et al1l2010l ). 


4 DISCUSSION 

We have studied the formation and future evolution of the 
first Be-BH binary system MWC 656. The study was carried 
out with population synthesis methods and it employed the 
standard model of single star and binary evolution. 

The formation of the system requires just two distinc¬ 
tive evolutionary steps. First, the massive primary (a BH 
progenitor) initiates a CE that is needed for the formation 
of close binary with an orbital period similar to the one ob¬ 
served for MWC 656. Then, a Wolf-Rayet star (an exposed 
core of the primary) explodes in type Ib/c SN and forms a 
~ 5M0 BH and makes the system eccentric, again as ob¬ 
served for MWC 656. At this point we note the formation of 
a massive binary consisting of a BH and a 10 — I6M0 B/Be 
star. 

Two competing scenarios are generally considered to 
explain the emission line phenomenon of Be stars. In both 
scenarios it is considered to be a direct result of rapid ro¬ 
tation and the presence of an outflowing disk. In the first 
scenario the B star is born as a substantially rapid rota¬ 
tor and retains this spin throughout its MS evolution. In 
the case of our evolutionary calculations it would also have 
to retain its spin throughout the CE phase initiated by its 
companion. In the second scenario the B star was spin up 
through interaction with the progenitor of the black hole, 
most likely through accreting a small amount of mass. Just 
a few percent of mass gain would be sufficient. 

In our model we do not observe significant accretion 
and spin up of the secondary star during the CE phase, 
giving some support to the first scenario to explain the Be 
phenomenon in this system. However, our rapid binary code 
relies on simplified prescriptions for stellar evolution and ap¬ 
proximate treatment of the stars in response to mass trans¬ 
fer. It would be worth further investigating these formation 
channels in more detail with a full binary evolutionary code 
which properly solves the structure equations for both stars. 
Particularly interesting promising possibilities for the for¬ 
mation of Be-BH systems include a phase of stable mass 
transfer (and thus spin-up by mass transfer) preceding the 
CE phases (delayed CE) or a phase of wind or atmospheric 
RLOF. 

At present the precise conditions for the Be phe¬ 
nomenon are not well understood. With our models we can¬ 
not make definite conclusions about the fraction of B stars 
that will show the phenomenon (and thus the fraction of 
B-BH binaries that appear as Be-BH binaries). In principle, 
explicit modeling of the stellar spins, the spin up process 
and the moments of i nertia will prov i de mo re i nformation, 
for ex ample as done bv ide Mink et all (l2013l) an d IShao fc Lil 
(l2014h under the assumption of rigid body rotation. How¬ 
ever, such simulations will not constrain the many remaining 
uncertainties. Among the most important unknowns is the 
mass transfer and the resulting spin up during accretion. It is 
still not well understood how the accretion stream interacts 
with the star or the accretion disk and how the accreting 
star responds when its outer layers are spun up. Another 





















The formation and evolution of MWC 656. 13 


important ingredient is the transfer of angular momentum 
throughout the star. Detailed stellar evolutionary models 
with and without interior angular mome ntum transport by 
magnetic fields give different results fe.g. iBrott et al.|l201ll : 
lEkstrom et al.ll2012l l. The situation is even less clear in the 
case of CE evolution. Particularly, it is unclear whether a 
star ca n accrete mass and spin up befo re or during the CE 
phase. (iMacLeod fc Ramirez-Rxiijl2015l l find that the inspi- 
raling object may increase its mass by a few percent. Thi s 
may be sufficient to spin up the accreting star dPacketlllQSlTl . 
As a general caveat the CE hydrodynamical simulations are 
still at the early st age and do not reproduce the observations 
dPassv et al.ll2012 '). Finally, the initial distribution of stellar 
rotation rates is not well constrained. The observations for 
young massive stars show a bimodal distribution with most 
stars rotating with slow or moderate rates with a tail of 


rapid rotat o rs reaching breakup velocity (e.g. Dufton_et_^ 


20n, I 2 OI. 1 I : lRami'rez-Agudei3 1201.11 : ISimon-Di'az fc Herrerc 


2014 . and references therein). Considering all these uncer¬ 


tainties we opted for not following the spin evolution explic¬ 
itly in this study. We investigated the formation channels 
of B-BH systems. A certain fraction of these will (intermit¬ 
tently) show the Be phenomenon. This fraction is uncertain. 
It may be as low as 0.3 if the ratio of Be stars to B stars is 
comparable to that observed in clusters. It may be higher if 
the interaction with the progenitor of the BH caused the Be 
phenomenon. A detailed investigation of the precise condi¬ 
tions for the occurrence of the Be phenomenon could prove 
to be very valuable. 

Although it was found that several tens of B-BH bina¬ 
ries are currently residing in the Galactic disk, we find that 
there is only small probability 1/100) of having system 
resembling MWC 656 in terms of component masses and or¬ 
bital period. In particular, in our simulations we find more 
systems with shorter orbital periods 5 — lOd), less mas¬ 
sive B stars 3 — 5M©) and slightly more massive BHs 
7 — 8Mq). For comparison, MWC 656 has orbital pe¬ 
riod of ~ 60d, Be star mass of 10 — 16M© and BH mass 
of ~ 5Mq. If MWC 656 is representative of the intrinsic 
Galactic population of Be-BH binaries, it means that ei¬ 
ther our employed evolutionary model of massive stars or 
our adopted approach to CE and/or BH formation needs 
to be re vised. We hav e empl oyed non-rotating stellar mod¬ 


els from iHurlev et ah J 2000|) with stellar winds corrected 
for clumping from Ivink et ahl (1200 il l. If a massive primary 
happens to be a fast rotator then our non-rotating models 
may not be a good choice, as one would expect smaller radii 
and larger cores for rapidly spinning stars. This would af¬ 
fect the development and the outcome of the CE phase and 
would al so most likely lead to the formation of mor e mas¬ 
sive BH (iLeitherer et al.l (l20l4 l. [de Mink et al.l (l2013ll '). The 
decreased stellar wind mass loss rates that we use are typi¬ 
cally adopted in most recent evolutionary studies of massive 
stars. However, it appears that there may be some observa¬ 
tional evid ence for higher wind mass loss rates from mas¬ 
sive stars jEldridge et al.l[2008l l. Had we adopted stronger 
winds, the primary core mass would decrease and to some 
extent it would counter-balance the evolutionary effects of 
fast rotation with the formation of a more massive BH. For 
C E evolution we have adopted the energy balance model 
of I WebbinkI Jl984j upda ted with the ph ysical estimates of 
primary binding energy dXu &: Lill2010bl lcJl. It appears that 


earlier claims that the core definition may change the post- 
CE binary separation by almost ~ 2 orders of magnitude 
(iDewi Si Taurid boOlf l for stars with M < 20Mq (NS pro¬ 
genitors) do es not apply to m ore massive stars (i.e., BH 
progenitors; IWong et al. (l2013lH . However, it is not at all 
clear that the energy balance model is good approximation 
for C E evolution, but no better model exists at the mo¬ 
ment dlvanova et ahllioisl l. For the B H formation we em - 
ploy the rapid supernova model from iFrver et al.l d2012|j . 
This model explains observe d mass gap between NSs and 
BHs dBelczvnski et al]l2012l l. the lack of compact objects 
in 2-5M© mass range. Our model also allows us to re- 
produce the Galactic a nd extragalactic BH mass spectrum 
dBelczvnski et al]l2010ll . It is possible that the mass gap is 
caused by some observati onal bias involved in BH mass mea¬ 
surements as proposed bv iKreidberg et ^ d2012|j . However, 
we have shown that our results do not depend sensitively on 
this aspect of evolution (i.e., BH formation mass; see model 
Vi). We have adopted an asymmetric mass ejection mecha¬ 
nism for BH natal kicks and this results in small BH kicks 
(decreased with the amount of fall back). The empirically 
derived kick velocities for 14 Galactic BH binaries are in¬ 
conclusive and the observa tional data allows for b oth small 
and high BH natal kicks dBelczvnski et alll2012| j. Had we 
adopted high BH natal kicks (similar t o those measured for 
single pulsars; e.g.. lHobbs et all d2005h l the formation rates 
of B-BH binaries would decrease by large factor (~ 1-2 or¬ 
ders of magnitude). In such the case we would expect very 
few Be-BH binaries to reside currently in the Galaxy. We 
have also tested the alternative model with BH natal kicks 
all set to zero. For this model we note a significant (a fac¬ 
tor of ~ 10; model V3) increase of MWC 656-like systems 
currently predicted to reside in Galaxy. 

With t he recent discovery of a Oe star with a compact 
companion dGlark et al.l l2015l ~) we can await more Be-BH 
binaries to be found. When more Be-BH binaries are identi¬ 
fied it will be essential to measure the intrinsic distribution 
of their orbital periods. Matching the intrinsic distribution 
with evolutionary models may allow us to improve our un¬ 
derstanding of BH natal kicks and may offer some insights 
into the inner workings of GE, as these processes seem to 
be the most important factors affecting the formation of 
MWG 656-like systems and B-BH binaries in general. 

The future evolution of MWG 656 offers a very 
interesting potential of forming a BH-NS system. 
In particular, close BH-NS systems are a class of 

J ravitational-radiation sources for Adv a nced LIGO 
Harry fc the LIGO Scientific Go llaborati^ I2OI0I I and 
Virgo ( Virgo GollaborationI I2OO9I 1 detectors, that are ex¬ 
pected to start operation in a few years. We have estimated 
the formation rate of close (delay time from formation on 
ZAMS to BH-NS coalescence shorter than 10 Gyr) BH-NS 
systems from binaries similar to that of MWG 656. We 
have translated the formation rate to detection rate of 
BH-NS mergers by advanced LIGO/Virgo. The advanced 
LIGO/Virgo detection rate is found to be up to 1 detection 
in 5 years. 

This empirically inferred detection rate is compara¬ 
ble to the rate obtained from analysis of Cyg X- 3 binary 
(1 detection in 10 years; iBelczvnski et al] (I 2 OI 3 I H and is 
much higher than obtained fo r Gyg X-1 (1 detection in 100 
years; IBelczvnski et al] (I 2 OIIIH . These empirical estimates 
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are based on existence of some particular binary stars and 
therefore they are only lower limits (i.e., formation of BH- 
NS only along one very specific formation channel) to the 
detection rate. However, it is encouraging that the observa¬ 
tional evidence, although indirect, increases to support the 
existence of close BH-NS systems. It is now the total of three 
systems: Cyg X-1, Cyg X-3 and MWC 656 that were shown 
to be potential BH-NS progenitors. Other known high-mass 
X-ray binaries w ere analyzed and exclu ded as potential BH- 
NS progenitors dBelczvnski et aHl2012l l. Recent population 
synthesis analysis of overall formation channels of double 
compact objects further supports the empirical evidence 
with estimates of BH-NS merger advanced LIGO/Virgo de¬ 
tection rat es at the level 0.03—5.7 yr“^ (from 1 in 30 years to 
6 per year; IPominik et al.l (l2014l'l'l. Broader det ection range 
is reported by Mennekens fc VanbeverenI ||2014| ~I: 0.04 — 484 
yr“^ (from 1 in 25 years to 484 per year). So far there are 
no known stellar-origin BH-NS binaries. At the moment the 
only known potential BH-NS binary consists of a central 
Galactic supermassive BH (Sagittarius A* ) and its nearby 
magnetar PSR J1745-2900 IIEatoughll2013l ). 


5 CONCLUSION 

At pr esent we know around 180 X-ray binaries dZiolkowskil 
l2014h . For ~ 120 of them the nature of compact object was 
confirmed to be a neutron star. Just in 2014 the first Be 
X-ray binary with a black hole was discovered (MWC 656; 
dCasares et al.|[2014h '). 

In this study we have investigated the possible evolu¬ 
tionary scenarios leading to the formation of B-BH systems, 
part of which are Be-BH binaries. It was found that the B- 
BH progenitors experience common envelope phase. For the 
majority (90%) of cases the CF donor is already an evolved 
star (CHeB), which increases chances of envelope ejection 
and a survival of the binary. In our simulations we used a 
code based on a non-rotating stellar models. Therefore, we 
are not able to resolve the issue of the origin of rapid rota¬ 
tion of Be stars. In our models it may be either connected 
to high initial star rotation or to spin up in CF (or pre-CF 
mass transfer) phase. 

We checked the sensitivity of our results on the various 
evolutionary parameters assumed in our simulations. For all 
presented models we expect up to few tens of B-BH systems 
to reside in the Galaxy at any given time, among which 
around 1/3 contain a Be star, but the chance to have a sys¬ 
tem very similar to the MWC 656 is less than a few percent. 
The future discoveries of Be-BH systems may allow us to de¬ 
termine the intrinsic distribution of orbital parameters and, 
by matching them with models, to improve our understand¬ 
ing of phases like CF and SN, which play an important role 
in the evolution of not only Be-BH binaries, but binaries in 
general. 

We investigated the fate of systems similar to 
MWC 656. The future evolution of such systems may lead to 
the formation of BH-NS binaries. We find that 18% (7%) of 
such population will form close (wide) BH-NS systems, while 
the rest merges in the second CF phase or is disrupted in 
the second SN. We have repeated the prediction for a sys¬ 
tem exactly like MWC656. Due to its favorable configura¬ 


tion, MWC656 is very likely (~ 77%) to form a close BH-NS 
system that will merger within 10 Gyr. 

These results make MWC 656 along with Cyg X-1 and 
Cyg X-3 the only reported potential progenitors of BH-NS 
binaries. The existence of MWC656 alone implies that the 
detection of BH-NS mergers by advanced LIGO/VIRGO can 
be as high as 1 in every 5 years. 
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APPENDIX A: CONSTRAINTS ON THE SN 
EXPLOSION IN MWC 656 


In this Appendix we report on the space velocity of MWC 
656 and on constraints on the SN explosion that originated 
the BH in this system. 

The position of MWC 656 in the ICRS (epoque J2000) 
according to Ivan LeeuwenI (l2007l l is: 

a = (22 42 57.30295 ±0.67) mas 

(Al) 

5 = (+44 43 18.2525 ±0.72) mas 


where a and 5 means right ascension and declination, re¬ 
spectively. Its proper motion from the same reference is: 

■ cos((5) = (—3.56 ± 0.72) masyr”^ 

_i (A2) 

+5 = (—4.05 ± 0.76) masyr 

From l^sares et al.l (|2014| ) the radial velocity of MWC 656 
is (—14.1 ± 2.1) kms~^ and its distance is (2.6 ± 0.6) kpc. 

The Galactic rotation curve of Model A5 in Table 4 of 
I Reid et ^ ll2014h provides: 


Rq = (8.34 ±0.16) kpc 
00 = (240 ±8) kms"^ (A3) 

deQ/diiQ =(—0.2 ±0.4) kms~^kpc~^ 

The same model provides a peculiar space velocity of the 
Sun relative to the Local Standard of Rest of: 

Uq = (10.7 ± 1.8) kms“^ 

R© = (15.6 ± 6.8) km s“^ (A4) 

Wq = (8.9 ± 0.9) kms~^ 


The coordinates adopted for the North Galactic Pole 
and Zero longitud e are those listed in the Appendix of 
iReid fc Brunthalerl il2004l l: 

apoie = 12'‘ 51™ 26.282“ = 192.8595083 deg 

= +27° 07 42.0l" = 27.12833611 deg (A5) 
e = 122.932 deg 


Using all these values we obtain a peculiar space velocity 
of MWC 656: 


1 ^spaceMWC 656 = (22.5 ± 14.9) ^m s ^ (A6) 

The uncertainty was calculated using values listed in 
Tab. which added in quadrature give 14.9 kms 
Therefore, the measurement of space velocity is only signifi¬ 
cant at the 1.5 sigma level. Even for a fixed distance with no 
uncertainty the total uncertainty is 12.7 kms“^, providing 
not even a 2-sigma detection. To obtain a 3-sigma detec¬ 
tion it is required to provide null uncertainties for both the 
distance and the galactic rotation curve. 

The masses and orbital parameters of MWC 656 from 
ICasares et al.l (|2014| ) are: 

Ml = (13±3)Mo 
Ma = (5.3 ± 1.5)M0 
P = (60.37 ± 0.04)d 
e = 0.10 ±0.04 


Table Al. The contribution of individual uncertainty to the 
space velocity uncertainty of MWC 656. 


Considered uncertainty 

Measured value (kms 

Proper motion 

5.2 

Radial velocity 

1.7 

Distance 

7.8 

Galactocentric distance 

0.5 

Galactic rotation of the Sun 

6.8 

Galactic rotation around MWC656 

7.0 

Peculiar space velocity of Sun 

5.8 


Table A2. The contribution of individual uncertainty to the un¬ 
certainty of the mass lost in the SN explosion of MWC 656. 


Considered uncertainty 

Measured value (Mq) 

Space velocity 

2.6 

Mass of Be star 

2.0 

Mass of BH 

0.6 

Re-circularized period 

0.02 


Considering the large orbital period and the small ec¬ 
centricity it is reasonable to assume that the current value 
of eccentricity is similar to the one just after the SN explo¬ 
sion, CpostSN. In such a case, the mass lost during the SN 
explosion is similar or slightly above to: 

Am = epostSN(Mi ■ M 2 ) = (1.8 ± 0.8)M© (A8) 

Using the formalism described in iNelemans et ahl lll999ll (al¬ 
though the orbital period here is much larger than 7 days) 
we find reduced mass, re-circularized period and initial or¬ 
bital period: 

M = 0.91 ±0.03 

Precirc = (59.5 ± 0.7)d (A9) 

Pinit = (49 ± 4)d 

With all these parameters we obtain an expected MWC 656 
space velocity of: 

t’expMwcese = (10.2 ± 4.7)kms"^ (AlO) 

to be compared to the measured space velocity of (22.5 ± 
14.9) hms”*^ from Eg. IA6I 

Both values are compatible at the 1 sigma level, indicat¬ 
ing that there is no need of any additional kick to produce 
the observed space velocity with a symmetric SN explosion 
with a mass loss of (1.8 ± 0. 8) M©. _ 

If we use the Eq.7 from lNelemans et al.1 lll999ll to com¬ 
pute the mass lost in the SN explosion considering the mea¬ 
sured space velocity and all other values with their corre¬ 
sponding uncertainties, we obtain: 

Am = (4.0 ± 3.4)M0 (All) 

where the uncertainty was calculated using values from Ta¬ 
ble IA2I As expected, this is compatible with the value de¬ 
rived from the current masses and the eccentricity of the 
orbit. 

In conclusion, everything seems compatible with no kick 
and a moderate mass loss of a few solar masses to pro¬ 
duce both the observed eccentricity and the space velocity 
of MWC 656. 




































